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Introduction 

The purpose of this research was to study the propagation of glactic ions through 
various materials. Galactic light ions result from the break up of heavy ion particles 
and their propagation through materials is modeled using the one-dimensional Boltzmann 
equation. When ions enter materials there can occur (i) the interaction of ions with orbital 
electrons which causes ionization within the material and (ii) ions collide with atoms 
causing production of secondary particles which penetrate deeper within the material. 
These processes are modeled by a continuum model. The basic idea is to place a conti ol 
volume within the material and examine the change in ion flux across this control volume. 
In this way on can derive the basic equations for the transport of light and heavy ions in 
matter. Green’s function perturbation methods can then be employed to solve the resulting 
equations using energy dependent nuclear cross sections. 
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Throughout this report we will use the following symbols and notation: 


— * 

<j>j(x, Q, E) is the flux of ions of type j moving in direction Q 

having units (#part.icles/cm 2 — sec — sr — Mev/amu) 


E 

Aj 

a j(E) 

Sj(E) 


R 3 (E) 

j 

Q 




n 

r 

X 

P 

Zj 


"i 


is the ion energy. (Mev/amu) 
is the atomic mass of the jth type ion (amu) 
is the macroscopic cross section, (cm -1 ) 
is the average energy loss per unit length or stopping power 
or linear energy transfer (Mev/cm). 

is the slowing down range for type j ions, (cm) Rj (E) = f Q s](E') 
is the ion type. 

is a unit vector in the direction of propagation. 

is production cross section of type j ions with energy E and direction Q 
by collision with type k ions of energy E' and direction ft 
having units of (cm — sr — Mev/amu) 
is the outward directed unit normal to boundary, 
is vector to boundary point, (cm) 

is the position vector to arbitrary point in region (cm) x = pEl + x n 

is the projection of x on Q (cm) 

is the component of x perpendicular to Q direction. 

charge of jth type ion. 

Z 2 

a parameter defined as Vj = jA . 


The basic Boltzmann equation results from examination of a control volume placed 
within the material. We find that 


/ Change in ion flux \ 
within a volume element 


V 


) 


( Gains within the'' 
volume element I — 


V 


/Losses due to any\ 
nuclear collisions 


\ 


This gives the Boltzmann equation 




=vf 


<P 0 (x,n,E) = / dE' j dVHa ik <f> k { x,ft,E r ) (1) 

k>j Je J 


The equation (1) is to be associated with the geometry of figure 1. 
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Figure 1. General Geometry for Boltzmann’s equation 
Multiply the equation (1) by Sj(E) and define the quantities 

4,(2 ,fi,E) = S j {E)(l> j (x,n,E) 

G, (x, ft, E) = Sj (E) dE ' J d ^' °jk4>k{x, E') 


k>j 


to obtain 


ft- V- 


SjjE) d 
Aj dE 


+ <7j(E) 


<t> j (x,n,E) = G j (x,n,E). 


( 2 ) 

( 3 ) 

( 4 ) 


Note that ft • ^ is the directional derivative in the direction ft and that 


dfy d<p 3 dRj _ d<i> 3 Aj 
~dE ~ dRj dE ~ dRj Sj{E ) 


so that the equation (4) can be written as 

d d , 

+ &j(E) 


4>j(x,n,E) = Gj(x,rt,E). 


( 5 ) 


[dp dRj 

Introduce the characteristic variables (17 j,£j) given by the transformation equations 

rjj = p — Rj(E ) £j = p + Rj{E ) (6) 
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where p = O • x. Also introduce the variables 

Xj (Vj ,£j) = 4>j + f n, #) 

g j (Vj ,£j)= Gj (pCl + x n , n, £) (7) 

Vj(Vj,Q = <^(£) 


/ 




Figure 2. Geometry for characteristic variables 


By the chain rule we have 

4-M + 4 and = + 

sp “ a% SR, 8% 1 ' 

so that the equation (5) simplifies to 

^2 7^- + cfj ^ Xj (r?j , £j) = 9j (Vj > ) 
in terms of the new variables. This equation can be integrated using the integrating factor 

( 9 ) 


(8) 


to obtain 


i r j — 

ex P g J °^^j) dr ] 

1 f Vj 

Xj in 3 ,&) = exp[- - / a, in' , Zj)drf]xj (a, £j) 

2Ja v ( 10 ) 
+ \J* exp[-i J aj{g'\$j)dri”}gj('n' ,£,j)dri’ 
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where a is any real number. Consequently, the solution to the equation (4) can be written 


as 


E) = f (a, p — Rj(E))<pj Q(£ i +a)fi + fn,a.R i 


rp-Rj(E) . i ^ c rff \ 

J a 


+ 


(ii) 


dr]' 1 


where 


f(a,p— Rj(E)) = exp 
From equations (6) we find that 


-u 


i rP—Rj(E)-_ c . L. 


( 12 ) 


2p = ?7 j + Cj and 2Rj{E) = £j - rjj 


(13) 


so that when rf = a we will have p = ^(a + £j). Observe from figure 2 that along the line 
of integration we will have = constant. The value of a is selected such that pQ + x n — T 
is a point on the boundary. Thus, the vector + x n = f dotted with (1 gives the 

value 

a = 2Q-?-Z j = 2d- p-Rj{E) (14) 

where d = • f . Note that when E = E' and rjj = r] 1 we have from equation (13) that 

2 Rj{E') = ^ - r)' (15) 


or 


e=R i'fi-r-) = R } 


1 f P + Rj( E ) — with dE > = ^£ 1 dr]' 


(16) 


and similarly by changing symbols when 


+ we have d E" = =W±dr]". 


e" = r: 


2A., 


2 j 

We examine the limits of integration in equation (11) and observe that when rj 
have 

2Rj{E') = p + Rj(E) — cl 

and from equation (14) we have 


(17) 


= a we 


(18) 


2d = p T Rj ( E ) T cl. 


(19) 
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Adding the equations (18) and (19) we find 


Rj (E 1 ) + d — p + Rj (■£") 


( 20 ) 


or 


E' = Rj{p -d + Rj(E)). 


( 21 ) 


Next we examine the lower limit of integration and find that wfien rf — p — Rj(E ), then 
2 Rj(E’) = p + Rj(E) - p + Rj(E) implies that E' = E. In the second term of equation 
(11) when rj" = a we again find that E" = RJ 1 (p + Rj{E ) — d ) and when rf' = p — Rj (E) 
then E u = E. Also, 

i(?j + n") = + {j - 2 Rj{E")) -ij- R A E ") = p+ r A e )- R A E '')- 

Consequently, the equation (11) can be written in the form 


<f>j (x, E) = Fj(E, RJ l (Rj{E) - d + p))^(r, 0, ^ {Rj(E) + p - d)) 


+ 


L 


RJ 1 (R j (E)+p-d ) 


Fj ( E , E")Gj{(p 4- A, (£) — Rj (E"))Q. 4- x n , f2, E ) ^ 




where 


Fj(El, E2) = exp 


_ t« A^V) 

Jei Si(E') 


Define the nuclear survival probability (reference Wilson 1977) as 


Pj (E) = exp 




A^ajE') , 
Sj(E') 


then the equation (23) can be written as 

FAE U E 2 ) = 


PjjE 2) 
P,(Ei)' 


Then from equation (22) we can write the solution to equation (3) in the form 

X f= * r\ rr\ _ Sj{Ej)Pj(Ej) ^ ^ ^ ^ ^ 

4 >j{x, n> e) - s ^ {E)p . {E) 4 >j{r,n,Ej) 


g 


3 CLE' WjlEL r dE " [dfte jk {E , t E , ')M 2 +WE)-R j &))fitf 
Sj(E)Pj{E) J E , J 


( 22 ) 


(23) 


(24) 


,E") 


(25) 


where Ej = R~ l (p + R 3 (E) - d ), x = x n + pQ and E' and E" have been interchanged. 
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In the one-dimensional straight ahead approximation Q is a unit vector in the direction 
of x with p = x, x n = 0,Tjj — x Rj(E),£j = x + Rj(E) and f = 0. (i.e. the origin 0 
moves to the boundary x — 0). The equation (25) then reduces to 


<f>j (x, E ) 


Sj(Ej)P i(Ej) 

Sj(E)Pj(E) 





„ A,P,(E') 

SjiEWE) 


rOO 

/ dE" a jk (E\ E")<t>k{x + Rj{E) - Rj(E'),E") 
J E f 


(26) 


where Ej is determined from x and E such that 


E, =R; l lx + B- 1 (Ej). 


(27) 


The solution given by equation (26) can be expressed in terms of Green’s function as 

poo 

<Pj(x,E) = T / G jk (x ) E,E 0 )M^E Q )dE 0 (28) 

k>j Jo 

where (f> k {0,E 0 ) = f k {E 0 ) are boundary conditions. Substituting the assumed solution 
given by equation (28) into equation (26) we obtain 

E, Eo)M0, Eo)dEo = J Q (| fa f(0 ’ Ej ’ E °^ 0 ’ Eo)dE ° 


+ 


yf 3 dE' 'If 'l, [ dE"<r jk (E',E")y f G ke {x+ R^-R^E'), E” , E o )M0> E 0 ) dE 0 . 

Je S j\ E ) P j\ E )J E f p Jo 


Note that when l = K we can equate like coefficients and find that Gj m (x, E, Eq) must 
satisfy the integral equation 


\ Sj( E j) P j( E j) ^ {C\ T? T? \ 

Gjm (x, E,Eo)- Gj m (0, Ej.Eo) 


+ V f El dE' r dE " a ^ E '' E ")°km^ + Rj(E) - E "> £ o) 

f Z'-JE s j( E ) p ji E ) Je ' 

fc > J 


(29) 


subject to the boundary condition Gj m (0, E, Eq) = 5j m S(E — Eq), where the value for Ej 
is determined from the inverse relation Ej = Rj l (x + Rj(E)). The Gj m terms are written 
using the Neumann expansion as a perturbation series 


G jm ( x, E, Eq) = £ G% (x, E , Eo) 


(30) 
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with leading term 


M°)( r p p ) — Sj(Ej)Pj{Ej) . §(E — Ef)) 

G jm (x,E,E 0 )- s .( E)p ^ E) tjmWj h o). 


(31) 


with Ej = Rj x (x + Rj(E)). Note that when x = 0 we have Ej = E so that G^(0, E, Eq) 
satisfies the above boundary condition. The higher order terms are determined from the 
recursive definition / 

G<” +1) (x, £,£*) = 

£ f E ’ dE' r dE"a jt (E', E")G™(x + R,(E) - Rj(E'), E", £„). (32) 

and must satisfy the boundary conditions G ( j^ 1 \x, E, Eq) = 0 for n = 0, 1, 2, .... In the 
special case n = 0 the equation (32) reduces to 

G$(x,E,E 0 ) = 

f Ej jtp* AjPj(E') f°° f, /p/ p//\ r P ^ 

E J E dE sMm L A ' o) 

where = x + i?j(E) - i?j(F) + fl fc (E"). (be. treat x + Rj{E) - Rj(E") as an 

x* value. See for example equation (27) .) Again we observe that when x = 0 we have 
Ej = E and so the boundary condition at x = 0 is satisfied. 

Cross Section assumption 1 

For interactions dominated by peripheral processes we use 


o jm (E',E")=o jm (E")t 5(E'-E'j 
so that the equation (33) becomes 

G™(x, E, Eo) = 

r E i 


(34) 


k 

where 


V'' f E} ji 7 / AjPj(E') f°° j /■ p'b tt 1 ' r^//\ Sk(E k )Pk(E k ) . SiEt — En) 

e/ e dE M ) ( ) s^mE ;r ) km { k o) 


E‘ k = R;\x + R,(E ) - Rj(E') + R k (E")). (36) 

We integrate with respect to E" and observe that the only nonzero term occurs when 
E" = E'. This gives 

<*>> *°> = E C dE ' ~ E ° )Sk - (37) 
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where 


Rk(E' k ) —x + Rj(E) — Rj{E ) 4- R k {E ) 
Rj(E r ) — Rk(E') =x + Rj(E) — R k {E k ). 


(38) 


We know that so that the above can be written as 

! 

(jt - l) Rk(E') —x + Rj(E) - R k {E’ k ) 
or (l - Ej R,(E’) =x + Rj(E) - R k (E' k ) 


Thus, we can write 


R,(E’) = . Vk A x + Rj(E) - R„(E' k )) 

\Vk - Vj\ 

or R k {E')= — -Ax + Rj(E) -Rk(E'k)). 

| Vk Vj | 

Differentiate the equation (39) with respect to E k to obtain 


(39) 


R' k (E')dE’ = r Xl—(-R’ k (E’ k ))dE' k or c IE' = - S -^ dE' k (40) 


| V k - Vj 

The equation (35) can then be written as 

rE' 


\ Uk U i' Sk(E') 


<«(-.'.*>, -E <41) 

k JE ki 

The only nonzero contribution comes when k = m and = Eq an d so equation (41) 
reduces to 


G<-|> (x, E, E 0 ) = 


hjm(x, E, Eo , E') if ^ ( Rm(Eo ) -x)< Rj(E ) < Rm(E 0 ) - x 


jm 


otherwise 


where 


and 


hj m (x 


FF E -i ... WW (g) 3 P "( £ °) 

, E, E o, £ ) s .( E ) Pj { E] \ Vm - Vj \ P m (E') 


E' = R7 1 


v. 


m 


3 v - Vj 


[x + R j (E)-R m (E o)] 


(42) 


(43) 


(44) 
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That is, when E' k = E 0 and k = m we have from the equation (38) that 

=x + R j (E)-R j (E')+R m {E') 

—x + Rj(E) - Rm{E 0 ) 

=x + Rj(E) — R m (Eo) 

= , — -A x +-Rj{E) - Rk(Eo)). 

Wm-Vj I 

Also from the transformation equations (13) the % , , Vj > variables are related through 

the range scale factors Vj and v k , where VjRj — v k Rk • This produces the relations 

Vk ~ £k = ~2R k = ~ ~ ~(Vj ~ £j)’ 

Vk V k 

Then from the equations 


Rm{E o) 
Rj(E') - Rm(E') 
1 - Rj(E’) 


Vy, 


Rj(E') 


+ Vj =€k + Vk = 2p 


Vj - Zj =7 r(% - &) 


V 


we find that by adding the equations (45) and (46) that 




Vj 


I 'k 


Vj 


and subtracting (46) from (45) we obtain 

2 ° = ( 1_ ^) %+ ( 1 + ?) ?t ' 

Interchanging j and k in the equations (47) and (48) we find that 


Vk 
tk = 


_ 3 


Vj + V k 


2v k 

Vk ~ Vj 
2v k 


Vj + 
Vj + 


v k - v 


2vk 

Vk T Vj 
2v k 




(45) 

(46) 


(47) 


(48) 


Then when rjj is a value rf lying between the constants — £,• and +£j, (See Figure 2(b)), 
we will have 

f Vj + Vk \ , , f V k — Vj ^ 


( •'j 

Vk = 1 — 


Zk = 


2v k 

Vk ~ Vj 
2v k 


M 

)”' + ( 


Vk - Vj 

2v k 

Vk + Vj 

2v k 


€j- 
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Changing k to m we find 


v m - Vj \ 


Sm — 


Note that the boundary condition G jm (0, E, Eq) = 5 jm 8(E - E 0 ) can be written in the 


Gj m (0, E, Eq) = 5j m 5(R :) Hfj) - Eq) = 5j m 8{^j - Rj(E 0 )) — 6(£ m R m (Eo)) 
so that when = R m (E 0 ) we have 






Using the equations (6) and (49) we now calculate the inequality which occurs in the 
equation (42) . From the equation (10) , with a = — £, we have the inequality — £j < 7]' < r]j 
which implies 


-£j <V' < Vj 


-x — Rj(E) < - Rm{E 0 )- 


V m Vj 


-X < 2l/m - R m (E 0 ) + Rj(E)- ( ) [x + Rj(E)) < X 

V m ~ Vj \ V m Vj J 


Vm ~t~ Vj 

Vm Vj 
Vm H" Vj 


{x + Rj{E)) <x - Rj(E) 


Vm 4 " Vj 
Vm Vj 


X — X < 


Vm Vj 


+ 1 - 


Vm, + Vj 
V m ~ Vj 


Rj(E) < X + 


Vm 4 " Vj 
V m ~ Vj 


VjX < u m Rm ( Eq ) UjRj(E) < Vm% 

—UjX > VjRj(E)—u m R m (Eo) > —v m x 
VmRm{E o) “ VjX > VjRj(E) > v m R m (Eo) v m x 

h R m (E a ) - x >R,(E ) > ^ (R m (E 0 ) - x) 

Vj Vj 

^ (R m (Eo) - x) <Rj(E) < ^R m {Eo) - X 


Cross Section assumption 2 

We start with equation (33) and assume o , j m {E' , E") has a Gaussian distribution of 


the form 


&jm(E ,E ) — <jj m (E ) r— — exp 

V 27T 


(E' - E" - e jm y 
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we can then write the equation (33) in the form 


pEj pOO 

G™(x t E,Eo) = dE' dE"F jk S k 

J E w 


( 50 ) 


where 


with 


AjPj{E^) . , w/x Sk{EjJPk{EjJ g \ 

Fjk S j (E)P j \E) 3 ^ ’ ^ S k (E")P k (E") i 


£(. = R k \x + Rj(E) — Rj(E') + Rk(E")) 


(51) 


(52) 


The integration of (50) is over the region illustrated in the figure 3 in the limit as T — * oo. 
In expanded form the equation (50) has the form 


G<&{x,E,E 0 ) = 
f E i . f°° .. AjPdE') S k (E' k )P k (E' k ) 

Je ^ Je' ^ S j {E)P ] {E) ajk S k {E")P k {E") 6km k 0) ' 


( 53 ) 


where 


& jk — &jk{E ) , nr— ex P 
k v "7T 


(£' - E" - ( jt )‘ 



Figure 3. Limits of integration for Green’s function term. 
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Note that for the first integration in the E" direction we have E' is a constant. 
Consequently, we let 


r = 


E" + e jk - E' 

— ± Tin 


with dr - 


dE' 


\/2A jk V^A jk 

The equation (53) can then be written in the form 

i 

GfJ l (x,E,E 0 ) = 
r E i r°° 

} E iE j.„ “simftOT 

- ^jk 


(54) 


AjPjiE') _ _ r 2 S k (E' k )P k {E' k ) 

c. 7^EE r 6kmS{E k ~ Eo >- 


(55) 


where r = y/2Aj k r — ej k + E' and 

E' k = R;\x + Rj(E) - Rj(E') - R k (f)). 


(56) 


This integral can be simplified by using one of the mean value theorems for integrals and 
written as 


G™(x, E, Eo) = 


AjPi(e') . ,sue t .W E D, 


(57) 


fe 

%/2A~" 




with r* = y/2A jk r* - e jk + E ' and 

= i^(x + Rj(E) - Rj(E') - i? fc (f,)) (58) 

where r* is some mean value in the interval ( ^^- , oo) and when = Eq, then E is a 
solution of the nonlinear equation 

R k (E 0 ) = x + i?j(£) - Rj(E') - R k (V2A jk r* - e jk + E') (59) 


provided E < E' < E y Consequently, we can write 


(i) 

jm 


1 

2 

0 


(x,E,E 0 ) = 
AlEAEA t 

s j (E)P j (E) a om 



ifE < E' < Ej 
otherwise 


(60) 


where E' is a solution of the nonlinear equation (59) , r* is some mean value and erfc is 
the complimentary error function. 
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Figure 4. New limits of integration for Green’s function term. 

Another viewpoint 

By interchanging the order of integration in equation (33) we obtain the limits of 
integration illustrated in the figure 4 and the equation (33) can be written as 

G^l{x,E,E 0 ) = [ 3 dE" [ dE'F jk 8 km + lim [ dE" f ' dE'F jk 5 km . (61) 

JE JE 1 - too JE j JE 

Observe that along the line £’ /, =coiistant . we have from equation (52) that 

dR k (E' k ) dE' k dRj(P) 

dE' k dE' dE' 

Ak dE r k _ _ Aj 
S k (E' k ) dE' ~ Sj(E') 

ip/ A m Sj(E) , 

" dE - - AMEl) dE “ 

Hence, when k — m and S km = 1, the equation (61) reduces to 


gV(x,E,E 0 ) = 


A m S j (E') 

AjSmiE'n) 


dE' + lim / dE" 


AmSj(E') , 
Flm A,S m (E- m ) ^ 
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(63) 


The limits in the above equation are determined as follows. Observe that when E' — E 
and k — m the equation (52) gives 

Rm{E' m ) =X + Rj(E) - Rj(E) + Rm(E ") 
or E'm = Rm( x + Rm{E")) = E m \ 
and when E' = Ej and k = m the equation (52) gives 

Rm(E'J =x + Rj(E) - RjiEj) + R m (E") 

. But RjiEj) = x + R 3 {E) so that E' m = E" . Also when E' = E" and k = m we obtain 
from the equation (52) that 

R m {E' m ) =x + Rj(E) - Rj(E") + R m {E") 

R m (E' m ) =x + Rj(E) - R,(E") + 

is m 

Yi 

Vm 


or 


Rm(E' m ) =X + Rj (E) + - lj RjiE") 

E' m = J R- 1 (x + Rj(E) + (2- ~ l) Rj(E")) = E m 3 


with 


E' = R~\x + Rj{E) - Rm(E' m ) + R m (E")). 

Using the properties of the Dirac delta function we find that the only nonzero contribution 
to the integral dE' m occurs when E' m = Eq. In this case the integral given by equation 
(62) simplifies to 




jm' 

i>E 


[“ dE " (E')Pj CgO , P m (Ep) _ f 

J E . dE Sj(E)PAE) Jmi ' ’ S m (E")P m (E") n 

f T ^„A m S j (E')P j (E’)_ , P„(Eq) 

rSSo J Ei d SjiEjPjlE) lm(E ’ > S m (E”)P m (E") h 


(64) 


+ 


where 


and 


E' = R~\x + Rj{E) - Rm(Eo) + R m (E ")) 

(65) 

Ho 

if£/ m 3 < Eq < E m i 
otherwise 

(66) 

Mo 

if E" <E 0 < E m i 

otherwise 

(67) 
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In the solution 




f E . r i,AmSj(E r )Pj(E) r" l 

J E , dE szmsr^ ■ 


Pm{E o) 


:/i 


( 68 ) 


' S m (E")Pm{E”Y 

f°° .^.AmSiiE^PjiE') ~ (w , _ Pm{Eo) __ f 

+ J Ei dE -sjmwr lmi ] s m (E'WE") } 

wc make the approximation that E' remains almost constant over the above intervals ot 
integration. We define the quantity 

1 E + Eq .... 

E[ = Rj l {x + rj{E) - Rm{E o) + Rmi^Ej + 2 

and use a mean value theorem for integrals to write the solution of equation (68) in the 
form 


G<P(x,E,Eo) = 


1 Am Sj{E\)PAE[)Pm{Eo) 

2 S j {E)P j {E)Sm{E* 1 )Pm{E 

1 AmSiWmEtiPmjE o) 
2S J (E)P J (E)Sm(E*)Pm(EZ 


+- 


Lf _L r e-P 

_ /*oo 

2 / _ r 2 

-Pt«L ‘ 


dr 


where 


Ej <E{ < E 

Ej <El < oo 
Ej + €jk ~ 

ri %/2Ajfc 

E + e jfc - 

\/ 2 Ajfc 


T2 = 


Using the error function 




and complimentary error function 


rOC 

erfc(x) = / e 

J X 


dr 


the above solution can be written in the form 


G#i(*,£,£b) = 

1 Am^j{Ej)Pj^ ^PmiEo) _ f ( er f(r 2 ) — erf(ri)) 

2 S j {E)P j (E)S m {E* 1 )Pm(Et) 


+ 


A m S 7 (£Q fi(£i)^m(£o) r erfc ( ri ) 
2 S j {E)P j {E)S m {E* 2 )Pm{E S) 
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